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ABSTRACT 

Aims. To obtain limits on the variation of the fine-structure constant a and the electron-to-proton mass ratio fi over 
different cosmological epochs. 

Methods. A new approach based on the comparison of redshifts of far infrared (FIR) fine-structure lines and low-lying 
rotational transitions in CO is proposed which is in principle more sensitive by a factor of ~10 compared to QSO metal 
absorption lines. 

Results. Estimations of the quotient F — a 2 / fi obtained for two distant quasars J1148+5251 (z — 6.42) and BR 1202- 
0725 = 4.69) provide AF/F = (0.1 ± 1.0) X 10" 4 and (1.4 ± 1.5) x 10" 4 , respectively. The obtained limits are 
consistent with no variation of physical constants at the level of ~0.01% over a period of 13 Gyr. 

Conclusions. Upcoming observations of quasars and distant galaxies in FIR fine-structure lines of different species and 
in CO low rotational lines at the SOFIA, HSO, and ALMA are expected to improve the current limit by, at least, an 
order of magnitude. 
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1. Introduction 

Possible time variation of coupling strengths and particle 
masses is now being discussed in regard to the accelerat- 
ing expansion of the Universe. Theoretical models impos- 
ing extra dimensions predict that dimensionless quantities 
like the fine-structure constant, a = e 2 /(Tic) « 1/137, the 

■ electron-to-proton mass ratio, fx = m c /m p , etc. depend 
on the scale length of extra dimensions in Kaluza-Klein 

, or superstring theories (for a review, see Garcfa-Berro et 

■ al. 2007). This scale factor may vary with cosmic time giv- 
ing rise to variations of fundamental constants which are 
defined in the combined 4D and extra-D space-time. Very 
different patterns from linear and slow-rolling to oscillating 
variations are considered in contemporary theoretical mod- 
els (e.g., Marciano 1984; Mota & Barrow 2004; Fujii 2005). 
But until now, the evidence for any changes in fundamental 
constants has not been unambiguously asserted. 

For example, laboratory measurements do not show 
variations of a at the level of a/a — (—2.6 ± 3.9) x 10~ 16 
yr _1 (Peik et al. 2006). A lower limit can be set by the fis- 
sion product analysis of a natural reactor in Oklo (in units 
of 10~ 17 yr" 1 ): -0.6 < a/a < 1.2 and -4 < a/a < 3 
reported by Gould et al. (2006) and Petrov et al. (2006), 
respectively, which approximately corresponds to the red- 
shift z ~0.14. At higher redshifts, molecular rotational and 
atomic resonance transitions in combination with observa- 
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tions of the Hi 21 cm hyperfine line were used to con- 
strain different combinations of physical constants (e.g., 
Varshalovich & Potekhin 1996; Drinkwater et al. 1998; 
Carilli et al. 2000; Murphy et al. 2001; Kanekar et al. 
2005; Tzanavaris et al. 2007). A new constrain on A/i//i 
at z = 0.68, A/i//i = —0.6 ± 1.9 ppm (la statistical error) 
was recently reported by Flambaum & Kozlov (2007) from 
the comparison of absorption-line spectra of NH3 with CO, 
HCO + , and HCN in radio frequency regiorQ. This value of 
A.fi/[i being linear extrapolated to z = corresponds to 
fi/fi = (—1 ± 3) x 10~ 16 yr _1 . Here we adopt a cosmology 
with H = 70 km s" 1 Mpc^ 1 , fl A = 0.7, and Ct m = 0.3 
(look-back time t z = 6.3 Gyr at z — 0.68). 

In optical spectra of QSOs, limits on Aa/a of compa- 
rable accuracy can be obtained from metal absorption line 
measurements. Special observations aimed at such measure- 
ments were performed at two different spectrographs - the 
HARPS at the ESO 3.6-m telescope (Chand et al. 2006) 
and the UVES at the ESO 8-m telescope (Levshakov et 
al. 2007a). The following values for Aa/a were obtained 
(Levshakov et al. 2007b): -0.12±1.79 ppm at z = 1.15 and 
5.66 ± 2.67 ppm at z = 1.84 (la errors including statistical 
and systematic parts are indicated) . Assuming linear varia- 
tion with cosmic time, one finds a/a = (—0.1 ±2.1) x 10 -16 
yr _1 and a/a — (5.5 ± 2.6) x 10~ 16 yr _1 , respectively. 

To test possible systematics in Aa/a measurements 
with the UVES spectrograph, we performed additional con- 
trol observations of asteroids (Molaro et al. 2007). Asteroids 
provide an accurate radial velocity reference at the level of 

1 Note that /j, — m p /m c in Flambaum & Kozlov (2007). 
1 ppm = 10" 6 . 
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1 m s -1 , or 0.03 ppm in units of Aa/cv . Since no systematic 
shifts larger than 1.5 ppm have been revealed, the positive 
signal at z — 1.84 may be real or induced by other yet 
unknown systematics. 

At higher redshifts, z > 3.5, the changes in fundamental 
constants can be probed through observations of distant 
galaxies and QSOs in far infrared (FIR) and radio ranges. 
Carbon, nitrogen and oxygen FIR fine-structure lines were 
suggested to probe the interstellar medium in galaxies at 
cosmological distances and to search for z > 10 objects 
(Petrosianet al. 1969; Loeb 1993; Stark 1997; Suginoharaet 
al. 1999; Blain et al. 2000; Boselli et al. 2002; Papadopoulos 
et al. 2004; Schaeres & Pello 2005; Nagamine et al. 2006). 

In this paper we report on the limits on a 2 / \i estimated 
for the first time from the [C n] 158 /im line combined with 
CO rotational lines detected in the spectra of the most 
distant quasar Jl 148+5251 at a redshift z — 6.42 (Maiolino 
et al. 2005) and in the northern component of a pair of 
AGN hosts BR 1202-0725 at z = 4.69 (Iono et al. 2006). 
The corresponding look-back time is t z = 12.9 Gyr and 12.5 
Gyr which is, respectively, 93% and 91% of the age of the 
Universe. 

2. Computational method 

In the nonrelativistic limit and for infinitely heavy point- 
like nucleus all atomic transition frequencies are propor- 
tional to the Rydberg constant, TZ « 109737.3 cm -1 . In 
this approximation, the ratio of any two atomic frequencies 
does not depend on any fundamental constant. Relativistic 
effects lead to the corrections to an atomic energy which 
are proportional to the product (aZ) 2 lZ, where Z is atomic 
number. Corrections accounting for the finite nuclear mass 
are proportional to filZ. For atoms these corrections are typ- 
ically much smaller than relativistic corrections, but they 
become important for molecules. 

To study the dependence of atomic frequencies on a it 
is convenient to expand them near the present-time value 
of a in the co-moving reference frame: 

uj z = uj + qx + . . . , x = (a z /a) 2 - 1 . (1) 

Here u> and oj z are the frequencies corresponding to the 
present value of a and to a change a — > a z presumably at 
epoch z. The parameter q (so-called g-factor) is individual 
for each atomic transition (Dzuba et al. 2002). 

If a is not a constant, the parameter x differs from zero. 
In linear approximation, |Aa/ev| <C 1, the corresponding 
frequency shift, Alu = lu z — ui, is given by: 

Auj/uj = 2QAa/u 1 (2) 

where Q = q/to is the dimensionless sensitivity coefficient, 
and Aq = a z — a. 

For distant objects such a frequency shift would cause 
an apparent change in the redshift 

Aw/w = -Az/(l + z), (3) 

where Az = z — z is the difference between an apparent red- 
shift z and its true value z. If tu' is the observed frequency 
from a distant object, then the true redshift is given by 

l + z = uj z /u)', (4) 

whereas the shifted (apparent) value is 

1 + z = oj/ui'. (5) 



Relativistic corrections grow with atomic number Z, but 
for optical and UV transitions in light elements they are 
small: Q ~ (&Z) 2 <C 1. For example, the UV Fen lines used 
in the Aa/a measurements have sensitivities |Q| ~ 0.03 
(Porsev et al. 2007). All other metal lines usually observed 
in quasar spectra have even lower values of Q except for 
some of the Fei or Znii transitions (Dzuba et al. 2002; 
Dzuba & Flambaum 2007), which are, however, too weak 
to be suitable for precise measurements of redshifts from 
quasar spectra. 

However, in the radio and FIR ranges one can signifi- 
cantly increase the sensitivity to a-variation by looking at 
transitions between the fine-structure levels. In the non- 
relativistic limit (a — > 0) such levels are exactly degen- 
erated. The corresponding transition frequencies ujf s are 
proportional to {aZ) 2 what means that for these transi- 
tions Q = 1, i.e. they are about 30 times more sensitive to 
changes in a than UV lines. This gives 

Awf s /w fe = 2Aa/a . (6) 

In order to utilize high sensitivity of the fine-structure 
transitions in differential measurements of a, we need an 
independent reference source. Good candidates for this pur- 
pose are low- lying rotational lines of CO: carbon monoxide 
is the second abundant molecule after H2, it is observed 
in many extragalactic objects^, and frequencies of the rota- 
tional transitions of CO (hundreds GHz) fall close to the fre- 
quency of the [Cll] fine-structure transition 2 P3/2 — >- 2 Pi/2 
at 1900.539 GHz. 

In a good approximation, the frequencies w ro t of the 
rotational lines of light molecules are independent of a, but 
sensitive to the electron to proton mass ratio /1: Lu Iot ~ [ilZ, 
i.e. 

Aw ro t/Wrot = A/i//Z . (7) 

If a and/or /i changes in course of cosmic time, the 
apparent redshifts for the fine-structure line(s) and for the 
rotational line(s) will differ. Using Eqs. ([3]), ([6]), and (JT]), 
one finds 

(Zrot ~ S is )/ (1 + z)= Az/(1 + z) = 2Aa/a - A^/fi , (8) 

where z is a reference redshift. 

By introducing the parameter F = a 2 / /i, we can rewrite 
Eq.© in the form 

Az/(1 + z) = AF/F . (9) 

One can also keep in mind that variation of fi is by far 
larger than that of a, as predicted by some GUT models 
(e.g. Langacker et al. 2002; Flambaum 2007). 

3. Results 

Now we apply Eq.© to the observations of [Cll] 158 /im 
and CO emission lines from the quasars Jl 148+5251 and 
BR 1202-0725 to obtain a limit on variation of the param- 
eter F. 

2 Emission lines from the rotational transitions of CO were 
detected in ~40 galaxies at z > 1 (for a review, see Solomon & 
Vanden Bout 2005), and recently at z = 5.0267 towards SDSS 
J033829. 31+002156.3 (Maiolino et al. 2007) and z = 5.7722 to- 
wards SDSS J0927+2001 (Carilli et al. 2007). 
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The spectral observations of the [Cn] 158 fim emis- 
sion line in the quasar J1148+5251 were carried out with 
the Institut de Radioastronomie Millimetrique (IRAM) 30- 
m telescope at the frequency of 256.1753 GHz covering a 
bandwidth of 1 GHz with 256 channels spaced by 4 MHz 
(Maiolino et al. 2005). At this frequency the resolving power 
of the telescope is 9.6 arcsec (Half Power Beam Width) and 
the 1 GHz bandwidth corresponds to 1170 km s _1 . The re- 
sulting spectral resolution and the noise in the coadded and 
rebinned spectrum were, respectively, 56 km s _1 and 0.3 
mK (2.8 mJy), leading to a 30% accuracy of the flux den- 
sity scale. The [C n] line was detected at a significance level 
of 8<7 for the total exposure time of 12.4 h. The redshift and 
the peak intensity of the [Cn] line are Zf s = 6.4189 ±0.0006 
and ii58 = 11.8 mJy. The reported error a z = 0.0006 corre- 
sponds to the uncertainty of the line position measurement 
of a v ~24 km s _1 , which is about one bin size in the [C n] 
spectrum at the Nyquist limit of 2 resolution elements. 

Observations of the CO (7— >6) and (6— >5) emission lines 
were obtained with the IRAM Plateau de Bure interferom- 
eter at the frequencies 108.724 GHz (the total integration 
time T cxp = 22 h) and 93.206 GHz (T cxp = 14 h), respec- 
tively (Bertoldi et al. 2003). At about 5 arcsec angular res- 
olution (5.7" x 4.1" at 3.2 mm) the CO emission line is un- 
resolved and coincides within the astrometric uncertainties 
of ±0.3 arcsec with the optical position of the quasar given 
by Fan et al. (2003). The coadded 3 mm data were rebinned 
to 64 km s" 1 (J = 7 -> 6) and 55 km s" 1 (J = 6 -> 5) 
resulting in the accuracy of the line position measurements 
of cr„ ~36 and 24 km s _1 , respectively. These uncertainties 
are again of a bin size in the reduced spectra. The redshifts 
and the peak intensities of the CO (7— >6) and (6— >5) lines 

are, respectively, zll~ &) = 6.4192 ± 0.0009, J (7 _ 6) = 2.14 

mJy, and z^ 5) = 6.4189 ± 0.0006, J (6 _ 5) = 2.45 mJy. 

Weighting the reported rotational redshifts with these 
peak intensities, one obtains the mean z m t — 6.4190 ± 
0.0005. We will take this value for the quasar's systemic 
redshift z. Using the reported redshift Zf s and the averaged 
Zrot, Eq.© yields AF/F = (0.1 ± 1.0) x 10" 4 . 

The second [C n] line was detected at Z[ s = 4.6908 to- 
wards the northern component of the quasar BR 1202-0725 
(Iono et al. 2006). The profile of this line is similar to the 
CO (5— >4) and (7— >6) lines seen at z rot = 4.6916 from the 
same component (Omont et al. 1999). 

The [C n] 158 /itm emission was observed with the 
Submillimeter Array interferometer (SMA, Ho et al. 2004). 
The total exposure time at a redshifted [C n] frequency of 
333.969 GHz was T cxp = 19.6 h, and the angular resolution 
was 3.4" x 2.7". The coadded spectrum was averaged using 
120 km s _1 bin size resulting in the rms noise of 7.5 mJy, 
or the signal-to- noise ratio S/N~ 3 (the peak flux density 
~23 mJy as shown in Fig. 1 in Iono et al.). Assuming the 
uncertainty of the [Cn] line position as ~l/4 bin size, one 
gets the error a v ~ 30 km s _1 . 

The CO J = 5 -> 4 line observed with the IRAM 
interferometer is detected at the ~ 5er confidence level 
(T exp ~ 16 h), and the CO J — 7 — * 6 line observed 
with the IRAM 30-m telescope - at the 3a confidence 
level. The angular resolution in the interferometric observa- 
tions was 5.0 x 2.5 arcsec, and velocity resolution of about 
60 km s _1 . The redshift of the northern component is 
Z5_4 = 4.6916. The uncertainty of this value is, probably, 
25-30 km s _1 , i.e. approximately one resolution element, 



considering rather noisy line profiles shown in Fig. 2 in 
Omont et al. The angular resolution of the 30-m telescope 
for the 2-mm beam is 17 arcsec. The error of the reported 
redshift zj^g — 4.6915 ± 0.001 corresponds to the radial 
velocity uncertainty of 53 km s _1 . 

Taking into account that the angular resolutions are 
similar in observations of the [C n] and CO J — 5 — » 4 
emission lines, we can use their redshifts Zf s = 4.6908 ± 
0.0006 and z rot = 4.6916 ± 0.0006 (both errors correspond 
to<7„ = 30 kms" 1 ) to calculate AF/F = (1.4±1.5) X 10" 4 . 

4. Discussion 

While comparing the redshifted frequencies of different 
species to measure hypothetical variations of physical con- 
stants, one must account for random Doppler shifts of the 
line positions caused by non-identical spatial distributions 
of species (referred to as the Doppler noise hereinafter) 
which can mimic non-zero signals in Aa/a or Afi/fi 
or in a combination of these quantities (e.g., Levshakov 
1994; Carilli et al. 2000; Bahcall et al. 2004). To quan- 
tify uncertainties induced by the Doppler noise a sample of 
(Aa/a) /(Afi/fi) measurements is to be collected. 

The main problem here is how to estimate the dis- 
persion of random velocity shifts a v for a given system 
of spectral lines. In case of a large sample size the value 
of <7 V can be found from the scatter of points. For a sin- 
gle measurement, a guess for a v comes from the compari- 
son with data on velocity differences between spectral lines 
of similar species in nearby clouds. Observations of local 
galaxies show that the intensity of [C n] is strongly corre- 
lated with the intensities of the low-lying rotational lines 
of CO, the fine-structure lines of [Ci] AA370, 609 /xm and 
[Oi] AA63, 146 fim (Malhotra et al. 2001), and the fine- 
structure line of [Nil] A205 /im (Petuchowski & Bennett 
1993; Abel 2006). However, the surface distribution of the 
[Cii] emission may not precisely follow the actual 12 CO 
contours (Stacey et al. 1985). The CO rotational lines, if 
optically thin, are emitted throughout the whole molecu- 
lar cloud. As for the [C n] emission, it is usually enhanced 
at the edges of the molecular cloud in the photodissocia- 
tion regions (PDRs). Additionally, diffuse gas from the Hn 
regions can contribute to the intensity of the 'PDR' lines 
(Kaufman et al. 1999). However, the impact from the dif- 
fuse gas decreases with increasing gas densities and drops 
from ~30% at n H ~ 1 cm" 3 to -10% at n H ~ 10 3 cm" 3 
(Kaufman et al. 2006). 

Because both considered high-redshift [C n] emitters 
belong to the giant molecular gas complexes with M — 
10 10 M o (Omont et al. 1999; Walter et al. 2004), the con- 
tribution to the integrated observed intensity of the [C n] 
emission from the H n regions is probably negligible. In par- 
ticular, Iono et al. (2006) note an excellent agreement be- 
tween the [Cii] and CO profiles at z = 4.69. At z = 6.42, 
Maiolino et al. (2005) favor nn ~ 10 5 cm' 3 to explain the 
observed luminosities in FIR, CO, and [Cii] lines. 

Our estimate of the value of the Doppler noise for the 
CO and [C n] emission lines at z — 6 is based on typical 
parameters of the interstellar molecular clouds with similar 
gas densities: n = 10 2 - 10 4 cm' 3 , M = 10 2 - 10 4 M Q , 
size L = 2 — 20 pc (Table 1 in Mac Low & Klessen 
2004) . Using the velocity dispersion and region size relation, 
cr„(km s~ 4 ) = 1.10L°' 38 (pc), empirically derived for molec- 
ular clouds by Larson (1981) and then studied in detail by 
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Solomon et al. (1987) and by other groups (e.g., Falgarone 
et al. 1992; Miesch & Bally 1994; Caselli & Myers 1995; 
Frieswijk et al. 2007), we find a v ^1 — 3 km s _1 . Thus, 
we may assume that if the [C n] and CO lines arise co- 
spatially, then the velocity offset between them does not 
exceed the velocity dispersion within the cloud which is of 
order of a few km s _1 . On the other hand, our derived 
limit \AF/F\ < 10 -4 corresponds to a v = 30 km s _1 which 
is the error of the line position measurement. From these 
values, wc arc confident that the quotient F — a 2 / fi re- 
mains constant at the level of 0.01% to a look-back time of 
13 Gyr. 

The reliability of the result obtained can be compared 
with previous radio observations of molecules at intermedi- 
ate redshifts. Kanekar et al. (2005) used observations of H I 
21 cm and OH 18 cm lines from two absorbers at z ~ 0.7 
and found A/// = (0.44 ± 0.36 sta t ± l-0 syst ) x 10" 5 . Here 
/ = g p (a / fx) , and g p is the proton gyromagnetic ratio. 
The systematic error was set from the assumption that the 
Doppler noise is 3 km s _1 . However, this seems to be too 
optimistic taking into account observations of the nearby 
diffuse/translucent gas clouds in our Galaxy. For instance, 
absorption profiles of the OH 18 cm and HCO + 3 mm lines 
studied at 140-240 m s _1 resolution towards extragalactic 
compact continuum sources (Liszt & Lucas 2000) exhibit, in 
general, similar complex structures over the radial velocity 
interval of ~10-20 km s _1 , while showing small differences 
in details on the scale of 1-2 km s _1 . However, comparison 
of HCO + absorption with H I 21 cm absorption towards the 
same targets reveals significant discrepancies. Figs. 6 and 7 
from Liszt & Lucas show that there are rather wide veloc- 
ity intervals (~20 km s" 1 ) where H I absorption is seen but 
HCO + is absent. It seems quite possible that the Doppler 
noise for a sample of H I 21 cm and molecular lines could be 
as large as 10 km s _1 . Another set of observations of local 
warm gas clouds within 15 pc around the Sun (Redfield & 
Linsky 2007) show that the distribution of cloud velocity 
difference has a width of ~10 km s _1 . This means that in 
lower resolution observations of quasar intervening systems 
the line-of-sight velocity differences of order 10 km s _1 may 
occur on a rather short linear scale of ^10-100 pc. 

The limit on A/i//i = —0.6 ± 1.9 ppm by Flambaum 
& Kozlov (2007) was estimated from the NH3 transitions 
near A = 1.26 cm and mm-wave lines of CO, HCO + , and 
HCN detected from one absorption-line system. In this case 
Aii/ ii — 0.289Aw/c, which means that the uncertainty of 2 
ppm corresponds to Av ~ 2 km s _1 . As discussed above, 
the systematic error from the Doppler noise can be of the 
same order or even larger. 

To avoid the influence of the Doppler noise on differ- 
ential measurements of physical constants, lines of only 
one element should be utilized. For instance, ground state 
FIR fine-structure transitions of [Ci] AA370,609 Lim, [Oi] 
AA63, 146 fim, [N 11] AA122, 205 fim, [O III] AA52, 88 fim, and 
[Fell] AA15, 26 /mi can be used to directly constrain Aa/a. 
The sensitivity of this procedure to a- variation will be dis- 
cussed in the following paper. 

Up to now, [C 1] was detected in five high- 2 objects 
where CO has been observed: H1413+117, J=l-0 line 
(Barvainis et al. 1997; Weifi et al. 2005) and J=2-l line 
(Weifi et al. 2003) at z = 2.557; IRAS F10214+4724 and 
SMM J14011+0252, J=l-0 line at z = 2.285 and 2.565, 
respectively (Weifi et al. 2005); APM 08279+5255, J=l-0 
line at z = 3.913 (Wagg et al. 2006); PSS 2322+1944, J=l- 



line at z = 4.120 (Pety et al. 2004). Unfortunately, low 
signal-to-noise ( < 10) and large uncertainties of the posi- 
tions of [Ci] and CO lines (^25 km s -1 ) prevent us from 
constraining relative changes in constants with an accuracy 
better than 10~ 4 . 

5. Conclusions 

In this letter we propose to use the FIR fine-structure lines 
of atoms and ions in conjunction with low-lying rotational 
lines of CO for probing the cosmic time evolution of the 
fundamental physical constants, namely the quotient F = 
a 2 jii. 

The reported constraints on AF/F are obtained for two 
[C 11] emitters recently discovered at z — 6.42 (Maiolino et 
al. 2005) and z = 4.69 (Iono et al. 2006). We found no 
evidence for the variability of F at the level of 10 -4 over 
a period of 13 Gyr. The statistical reliability of this limit 
seems to be high enough since the expected value of the 
Doppler noise for the co-spatially distributed [C 11] and CO 
emission lines, a v ~ a few km s _1 , is less that the error of 
the line position measurement, a v — 30 km s _1 . However, 
to be completely confident it should be tested whether the 
empirical relationship between the velocity dispersion and 
region size derived on base of local measurements is also 
applicable to high-z galaxies. 

The result obtained is to be compared with other tests 
at extremely high redshifts: bounds from the height and 
the positions of the first and subsequent acoustic peaks of 
the spectrum of the cosmic microwave background radi- 
ation fluctuations [z ~ 1500, look-back time 13.8 Gyr), 
and from the relative abundances of light elements pre- 
dicted in the big bang nucleosynthesis (z ~ 10 10 ) restrict 
the fine-structure constant variations at the level of only 
2-3% (Ichikawa et al. 2006; Dent et al. 2007; Mosquera et 
al. 2007). Thus, the proposed method ensures a significant 
gain in the accuracy and works at the time scales compa- 
rable with the total age of the Universe. 

The limits on AF/F can be improved by measuring the 
line positions with better precision - as high as 1-3 km s" 1 
to be comparable with the expected Doppler noise. The 
advantage of radio observations for such precise spectral 
measurements is twofold: firstly, distant objects inaccessi- 
ble to optical observations can be probed, and, secondly, 
spectral resolution available in cm- and mm-range is much 
higher than in the optical. Additionally, the sensitivity to 
a-variation in FIR band is further increased by roughly a 
factor of 30 due to larger Q-factors, as compared to opti- 
cal band. In conclusion, the method proposed in this letter 
holds the promise of a higher accuracy than obtainable at 
present with QSO metal absorption lines. 

It is to note, that observations of both [C11] and CO 
lines in quasars and high-redshift galaxies will be one 
of the key tasks at the the Stratospheric Observatory 
for Infrared Astronomy (SOFIAj^, the Herschel Space 
Observatory (HSOj3, and the Atacama Large Millimeter 
Array (ALMA0 which is widely discussed in the litera- 
ture. Thus, we can expect that these upcoming observa- 
tions will contribute much to clarifying whether fundamen- 
tal constants vary with cosmic time or not. 

http://www.sofia.usra.edu/ 
4 http://herschel.esac.esa.int/ 
http://www.eso.org/alma/ 
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